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Abstract. Using a generalized parton model approach including spin and intrinsic parton 
motion effects, and assuming the validity of factorization for large pr jet production in hadronic 
collisions, we study the azimuthal distribution around the jet axis of leading pions, produced in 
the jet fragmentation process. We identify the observable leading-twist azimuthal asymmetries 
for the unpolarized and single-polarized case related to both quark and gluon-originated jets. 
We account for all physically allowed combinations of the transverse momentum dependent 
(TMD) parton distribution and fragmentation functions, with special attention to the Sivers, 
Boer-Mulders, and transversity quark distributions, and to the Collins fragmentation function 
for quarks (and to the analogous functions for gluon partons). 



1. Introduction 

Transverse single-spin and azimuthal asymmetries in high-energy hadronic reactions have raised 
a lot of interest in the last years (see e.g. Refs. [HE] and references therein). In particular, the 
huge spin asymmetries measured in the inclusive forward production of pions in high-energy 
pp collisions, at moderately large transverse momentum, cannot be explained in the realm 
of leading- twist (LT) perturbative QCD (pQCD), based on the usual collinear factorization 
theorems. 

Out of the theoretical approaches proposed in order to account for these measurements, in 
the following we will adopt the so-called transverse momentum dependent (TMD) formalism, 
which takes into account spin and intrinsic parton motion effects assuming a pQCD 
factorization scheme. Single-spin and azimuthal asymmetries are generated by TMD polarized 
partonic distribution and fragmentation functions, among which the most relevant from a 
phenomenological point of view are the Sivers distribution O |3] and, for transversely polarized 
quarks, the Boer-Mulders distribution [5J and the Collins fragmentation function [6j (similar 
functions can be defined for linearly polarized gluons, see e.g. Ref. [7]). 

Along the lines of [5], we consider the process p^^'^p jet + ir + X, presently under active 
investigation at RHIC, where one observes a large pT jet and looks for the azimuthal distribution 
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of leading pions inside the jet. A very preliminary version of this study was first presented in 
Ref. [9]. A similar analysis was performed in Ref. [TO], which however considered intrinsic parton 
motion only in the fragmentation process, drastically reducing the possible contributions to the 
asymmetry. Indeed, in that case, only the Collins effect for quarks is at work. In fact, Ref. [lOj 
aimed at studying only the Collins fragmentation function (FF), which should be universal, in 
a more simplified theoretical scheme for which factorization has been proven. Our approach 
is different in some respects. It is more general and has in principle a richer structure in the 
observable azimuthal asymmetries, since intrinsic motion is taken into account in the initial 
hadrons also. However, since factorization has not been proven in this case, but is rather taken 
as a reasonable phenomenological assumption, the validity of the scheme and the universality 
of the TMD distributions involved require an even more severe scrutiny by comparison with 
experimental results. On the other hand, at the present theoretical and experimental stage, 
we believe that combined phenomenological tests of different approaches are required to clarify 
the validity of factorization and, related to this, the relevance of possible universality-breaking 
terms for the TMD distributions. 

The plan of this contribution is as follows. In Sec. [2] we will briefly summarize the TMD 
generalized parton model approach and give the expression of the polarized cross section for the 
process of interest. In Sec. [3] we will present phenomenological results for azimuthal asymmetries 
discussed in the kinematical configuration of the RHIC experiments. Sec. 4 contains final remarks 
and conclusions. 



2. Formalism 

We denote with A and B two spin 1/2 hadrons (typically, two protons), with hadron B 
unpolarized and hadron ^ in a pure transverse spin state described by the four- vector Sa- 
Within a generalized TMD parton model approach, the invariant differential cross section for 
the process A{Sa) B ^ jet + tt + X can be written, at LT in the soft TMD functions, as 

Ejda /■ dXadXfo 2, j2i a/A,SA r t 1 \ 

d3p,d.d^fc^. = E /l^^^2^dfc^.dfc^,p,^,, AM,5.(Xa,fc^.) 

where and pj are respectively the energy and three- momentum of the observed jet. We sum 
over all allowed partonic processes contributing to the physical process observed, and {A} stays 
for a sum over all partonic helicities. Xa^b and k±a,b are respectively the initial parton light- 
cone momentum fractions and intrinsic transverse momenta. Analogously, z and k±T^ are the 
light-cone momentum fraction and the transverse momentum of the observed pion inside the jet 
w.r.t. the jet (parton c) direction of motion. 

All information on the polarization state of the initial parton a is contained in 
p'x'x'^'^ fa./A,SAi^a.^^-i-a)^ ^j^-^j^ depends in turn on the (experimentally fixed) parent hadron 
A polarization state and on the soft, nonperturbative dynamics encoded in the eight leading- 
twist polarized and transverse momentum dependent parton distribution functions. p'\'\f'^ 
is the helicity density matrix of parton a. Analogously, the polarization state of parton b 
inside the unpolarized hadron B is encoded into pi^^/ fb/Bi^b, k±b). The Ma a a a 's are the 

1} f) ' I) 

pQCD leading-order (LO) helicity scattering amplitudes for the hard partonic process ab — t- cd. 
The DJ^ A' (-^' ^^7r)'s are the soft leading-twist TMD fragmentation functions describing the 
fragmentation process of the scattered (polarized) parton c into the final leading pion inside the 
jet. More details can be found in Ref. |8J. 

We work in the AB hadronic cm. frame, with hadron A moving along the +Zcm direction, 
and define {XZ)cm as the production plane containing the colliding beams and the observed jet. 



with (pj)xcm > 0- t'^is frame Sa = (0, cos , sin (/)5^ , 0) and pj = r (cosh r/j, 1,0, sinhr/j), 
where rjj = — log[tan(0j/2)] is the jet (pseudo)rapidity. 

The calculation is performed by summing explicitly over all helicity indexes and inserting the 
appropriate expressions for the helicity density matrices of partons a, b and for the polarized 
distribution and fragmentation functions. After factorizing explicitly all azimuthal dependences, 
including those coming from the hard-scattering helicity amplitudes, collecting them and using 
symmetry properties under k±a^b — ^ —f^±a,b [8], one gets the final expression for the single 
transverse polarized cross section. This will have the following general structure: 

2da{(j)s^ , 0f ) ~ do-Q + dAcJo sin + dcJi cos (p^ + da2 cos 20^ + dAaf sin((/)5^ - (j)^ ) 

+ dAa+ sm{(Ps^ + dAa^ M^Sa " 2</'f ) + dAa+ sin(0s^ + 20^) , (2) 

where (f)^ is the azimuthal angle of the pion three-momentum around the jet direction of motion, 
as measured in the fragmenting parton helicity frame. The latter frame is related to the hadronic 
cm. frame by a simple rotation by 6j around Ycm = Vj [8\. 

In terms of the polarized cross section in Eq. ([2]), we can define average values of appropriate 
circular functions of (psA 'Pn ^ ™ order to single out the different contributions of interest: 

( WicPsA , ) ) (P„ fcx.) ;d^.,d0^d.fe,0.^) • 

Alternatively, for the single spin asymmetry we can, in close analogy with the case of semi- 
inclusive deeply inelastic scattering (SIDIS), define appropriate azimuthal moments, 

J dcPs^dcP^ [dai^sA^cPi^) + daicPsA + vr, <A^)] ' ^ ' 

where W{cpSj^, cp^) is again some appropriate circular function of cpSj^ and cp^ . 



3. Phenomenology 

In this section we present and discuss some phenomenological implications of our approach for 
the unpolarized and single-transverse polarized cases in kinematical configurations accessible 
at RHIC by the STAR and PHENIX experiments. We consider both central (?7j = 0) and 
forward (r/j = 3.3) (pseudo)rapidity configurations at a cm. energy ^/s = 200 GeV (different 
cm. energies, namely ^/s = 62.4 and 500 GeV, are also studied in [8j), aiming at a check of 
the potentiality of the approach in disentangling among different quark and gluon originating 
effects. 

We will first consider, for vr^ production only, a scenario in which the effects of all TMD 
functions are over-maximized. By this we mean that all TMD functions are maximized in size 
by imposing natural positivity bounds (and the Soffer bound for transversity |lHll2j): moreover, 
the relative signs of all active partonic contributions are chosen so that they sum up additively. 
In this way we set an upper bound on the absolute value of any of the effects playing a potential 
role in the azimuthal asymmetries. Therefore, all effects that are negligible or even marginal in 
this scenario may be directly discarded in subsequent refined phenomenological analyses. 

As a second step in our study we consider, for both neutral and charged pions, only the 
surviving effects, involving TMD functions for which parameterizations are available from 
independent fits to other spin and azimuthal asymmetries data in SIDIS, Drell-Yan, and e~^e~ 
processes. 



For numerical calculations all TMD distribution and fragmentation functions will be taken 
in the simplified form where the functional dependences on the parton light-cone momentum 
fraction and on transverse motion are completely factorized, assuming a Gaussian-like flavour- 
independent shape for the transverse momentum component. Concerning the parameterizations 
of the transversity and quark Sivers distributions, and of the Collins functions, we will consider 
two sets: SIDIS 1 [SI E] and SIDIS 2 [El [16]. Notice that the almost unknown gluon 
Sivers function was tentatively taken positive and saturated to an updated version of the 
bound obtained in Ref. [T7| by considering PHENIX data for the tt" transverse SSA at mid- 
rapidity production in polarized pp collisions at RHIC |18j . Furthermore, for the usual collinear 
distributions, we adopt the LO unpolarized set GRV98 and (for the Soffer bound) the 
corresponding longitudinally polarized set GRSV2000 [20j. For fragmentation functions, we will 
adopt two well-known LO sets among those available in the literature, the set by Kretzer |21j 
and the DSS one [22]. Our choice is dictated by the subsequent use of the two available 
parametrization sets for the Sivers and Collins functions in our scheme, that have been extracted 
in the past years by adopting these sets of FFs. 

Since the range of the jet transverse momentum (the hard scale) covered is significant, we take 
into account proper evolution with scale. Concerning transversity, in the maximized scenario 
we will fix it at the initial scale by saturating the Soffer bound and then letting it evolve. On 
the other hand, the transverse momentum component of all TMD functions is kept fixed with 
no evolution with scale. Notice that at this stage evolution properties of the full TMD functions 
are not known. 

In all cases considered, since we are interested in azimuthal asymmetries for leading particles 
inside the jet, we will present results obtained integrating the light-cone momentum fraction of 
the observed hadron, z, in the range z > 0.3. 

3.1. Azimuthal asymmetries in pp — )■ jet + it + X 

The symmetric part duo in Eq. ([2]) gets contributions by the usual unpolarized term, already 
present in the collinear approach, and by an additional term involving a Boer-Mulders (E) Boer- 
Mulders convolution for the initial quarks (or the analogous terms involving linearly polarized 
gluons); however even in the maximized scenario this last contribution is always negligible in all 
the kinematical configurations considered, hence we will not discuss it anymore in the sequel. 

In Fig. [1] we show the maximized (cos (j)^) and (cos 20^) asymmetries for vr"^ production in the 
central (left panel) and forward (right panel) rapidity regions as a function of pj-r, from pj^ = 2 
GeV up to the maximum allowed value, adopting the Kretzer FF set. Similar results are obtained 
using the DSS set. The coscj)^ asymmetry is generated by the quark Boer-Mulders (S" Collins 
convolution term, involving a transversely polarized quark and an unpolarized hadron both in 
the initial state and in the fragmentation process. The cos 2(f>^ asymmetry is related to the 
term involving linearly polarized gluons and unpolarized hadrons both in the initial state and 
in the fragmentation process, that is the convolution of a Boer-Mulders-like gluon distribution 
with a Collins- like gluon FF. Even the maximized contribution is practically negligible in the 
kinematical configurations considered. 

3.2. The bivers asymmetry Aj^ p'p — )• jet + vr + A 

In Fig. [2]we show the total observable Sivers asymmetry, and the corresponding quark and gluon 
contributions for vr"^ production, in the maximized scenario and adopting the Kretzer FF set, 
as a function of p^x in the central (left panel) and forward (right panel) rapidity regions. The 
maximized potential Sivers asymmetry can be very large in both cases. In the central rapidity 
region, the asymmetry is dominated by the gluon contribution at the lowest pj t range while gets 
comparable quark and gluon contributions in the large p^T range. A large Sivers asymmetry 
around p^T = 4-^6 GeV could then be a clear indication for a sizable gluon contribution. In 
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Figure 1. Maximized quark-originated {coscj)^) and gluon-originated {cos2cl)^) asymmetries 
for the unpolarized pp — t- jet + vr^ + X process at y/s = 200 GeV in two different rapidity regions, 
adopting the Kretzer FF set. 



the forward rapidity region, on the contrary, the quark and gluon contributions are comparable 
at low pjT values, while the maximized asymmetry is dominated by the quark contribution for 
> 4 GeV. Therefore, a large Sivers asymmetry in this kinematical range could be ascribed 
unambiguosly to the quark Sivers effect. 

In Fig. [3] we show, for both neutral and charged pions, the quark and gluon contributions to 
the Sivers asymmetry, obtained adopting the parametrization sets SIDIS 1 and SIDIS 2, and the 
updated version of the bound on the gluon Sivers asymmetry derived in Ref. |17j . in the forward 
rapidity region, as a function of pjT- The dotted black vertical line delimits the region xp ^ 0.3, 
with xp = 2pjL/\^, beyond which the SIDIS parameterizations for the quark Sivers distribution 
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Figure 2. Maximized total, quark-originated and gluon-originated Sivers asymmetries for the 
p^p — )■ jet -|- 7r+ -|- X process, at y/s = 200 GeV in two different rapidity regions, adopting the 
Kretzer FF set. 
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Figure 3. The estimated quark and gluon contributions to the Sivers asymmetry for the 
•p^p — )• jet + 7r + X process, obtained adopting the parametrization sets SIDIS 1 and SIDIS 2, at 
forward rapidity and \fs = 200 GeV. The dotted black vertical line delimits the region xp 0.3. 



are extrapolated outside the x region covered by SIDIS data and are therefore plagued by large 
uncertainties. This reflects on the fact that below this limit the two sets give comparable results, 
while above it they differ remarkably. Therefore, a measurement of this asymmetry might help 
in clarifying the behaviour of the quark Sivers distribution in the large x region, which plays 
a fundamental role for forward pion production at RHIC, and is not covered by present SIDIS 
data from HERMES and COMPASS experiments. 

3.3. The Collins (-like) ("^s^T^^^ ) ("^s^T 4'^ ) j asymmetries in p^p — t- jet + vr + X 

The quark generated asymmetry j^^^^'^^a^'I'^ ) comes from two distinct contributions: one 
involving the convolution between the term of the TMD transversity distribution suppressed 
in the collinear configuration and the Collins function; another term involving the convolution 
of the Sivers and Boer-Mulders distributions for the initial quarks with the Collins function 

for the final quark [an analogous term appears also in the Aj^ ^'^^^ '^^ ^ asymmetry]. We have 
explicitly checked that for the kinematical configurations under study both these contributions 
are always negligible already in the maximized scenario. Therefore we will not consider the 
sin((/>s^ + (j)^ ) asymmetry in the sequel. A similar situation holds also for the gluon generated 

^7V '^'^^^'^ '^^ asymmetry, where two contributions analogous to the quark ones discussed above 
but for linearly polarized gluons are involved. 

In Fig. m we present the quark ^'^^a 't'^ ) Collins asymmetry and the gluon A^^ ^'^^a '^^ ) 
Collins-like asymmetry in the maximized scenario in the central (left panel) and forward (right 
panel) rapidity region as a function of pjT, from pj-p = 2 GeV up to the maximum allowed value, 
adopting the Kretzer FF set. In the central rapidity region the quark Collins asymmetry is very 
small at the lowest pjT values, then increases almost linearly reaching about 8% at the upper 
range. Instead, in the forward rapidity region the asymmetry is (potentially) always large and 
increases almost linearly from about 25% to about 70% going from the lowest to the largest t 
values. Concerning the gluon Collins-like asymmetry, both in the central and in the forward 
rapidity regions it is of the order of 5% at the lowest t values, then starts decreasing slowly 
and becomes negligible at large pjT values. Similar results hold when adopting the DSS set. 

We consider now, for both neutral and charged pions, numerical results for the quark Collins 
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Figure 4. Maximized quark and gluon Collins (-like) asymmetries for the p^p — t- jet + vr^ + X 
process, at ^/s = 200 GeV in two different rapidity regions adopting the Kretzer FF set. 



asymmetry obtained adopting the parameterizations SIDIS 1 and SIDIS 2 for the transversity 
distribution and the Collins fragmentation function (no parameterizations are available yet in 
the analogous gluon case). It turns out that in the central rapidity region all the estimated 
asymmetries are practically negligible. Concerning the forward rapidity region, our results are 
shown in Fig. [SJ The Collins asymmetry for neutral pions results to be almost vanishing. 
For charged pions, similarly to the case of the Sivers asymmetry, the two parameterizations 
give comparable results only in the t region where the transversity distribution is reasonably 
constrained by SIDIS data (see the dotted black vertical line). A measurement of this asymmetry 
would be then very important and helpful in clarifying the large x behaviour of the quark 
transversity distribution. 
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Figure 5. The estimated quark Collins asymmetry for the p^p — t- jet + tt + X process, obtained 
adopting the parameterizations SIDIS 1 and SIDIS 2 respectively, at ^/s = 200 GeV in the 
forward rapidity region. The dotted black vertical line delimits the region xp ~ 0.3. 



4. Conclusions 

We have presented a study of the azimuthal asymmetries measurable in the distribution of 
leading pions inside a large-pr jet produced in unpolarized and single-transverse polarized proton 
proton collisions for kinematical configurations accessible at RHIC. To this end, we have adopted 
a generalized TMD parton model approach with inclusion of spin and intrinsic parton motion 
effects both in the distribution and in the fragmentation sectors. 

In contrast to inclusive pion production, where the Sivers and Collins mechanisms cannot be 
separated [S], and in close analogy with the SIDIS case, the leading-twist azimuthal asymmetries 
discussed above allow one to discriminate among different effects by taking suitable moments 
of the asymmetries. In principle, quark and gluon originating jets can also be distinguished, 
at least in some kinematical regimes. Hence, the proposed phenomenological analysis could be 
very helpful, for example, in clarifying the role played by the quark(gluon) Sivers distribution 
and by the Collins (-like) fragmentation function in the sizable single spin asymmetries observed 
at RHIC for forward pion production. At the same time it will give us the opportunity of testing 
the factorization and universality assumptions, and of gaining information on the size and sign 
of the TMD functions discussed, in a kinematic region not covered by SIDIS data. 

We finally stress that the unambiguous measurement of any of the asymmetries, other than 
the Collins one, discussed above would be a clear indication of the role played by intrinsic parton 
motion in the initial hadrons for the spin asymmetry sector in polarized hadronic collisions. 
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